home *** CD-ROM | disk | FTP | other *** search
-
-
-
- ____CCCCBBBBAAAALLLL((((3333FFFF)))) ____CCCCBBBBAAAALLLL((((3333FFFF))))
-
-
-
- NNNNAAAAMMMMEEEE
- CBAL, SCBAL - EISPACK routine. This subroutine is a complex version
- of BALANCE.
-
-
- SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ccccbbbbaaaallll((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh
- ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,, aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
-
- ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssccccbbbbaaaallll((((nnnnmmmm,,,, nnnn,,,, aaaarrrr,,,, aaaaiiii,,,, lllloooowwww,,,, iiiigggghhhh,,,, ssssccccaaaalllleeee))))
- iiiinnnntttteeeeggggeeeerrrr nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh
- rrrreeeeaaaallll aaaarrrr((((nnnnmmmm,,,,nnnn)))),,,, aaaaiiii((((nnnnmmmm,,,,nnnn)))),,,, ssssccccaaaalllleeee((((nnnn))))
-
-
-
- DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
- This subroutine balances a COMPLEX matrix and isolates eigenvalues
- whenever possible.
-
- On INPUT
-
- NNNNMMMM must be set to the row dimension of two-dimensional array parameters
- as declared in the calling program dimension statement.
-
- NNNN is the order of the matrix.
-
- AAAARRRR and AI contain the real and imaginary parts, respectively, of the
- complex matrix to be balanced. On OUTPUT
-
- AAAARRRR and AI contain the real and imaginary parts, respectively, of the
- balanced matrix.
-
- LLLLOOOOWWWW and IGH are two integers such that AR(I,J) and AI(I,J) are equal to
- zero if
- (1) I is greater than J and
- (2) J=1,...,LOW-1 or I=IGH+1,...,N.
-
- SSSSCCCCAAAALLLLEEEE contains information determining the
- permutations and scaling factors used. Suppose that the principal
- submatrix in rows LOW through IGH has been balanced, that P(J) denotes
- the index interchanged with J during the permutation step, and that the
- elements of the diagonal matrix used are denoted by D(I,J). Then
-
- SSSSCCCCAAAALLLLEEEE (J) = P(J), for J = 1,...,LOW-1
- = D(J,J) J = LOW,...,IGH
- = P(J) J = IGH+1,...,N. The order in which the interchanges
- are made is N to IGH+1, then 1 to LOW-1. Note that 1 is returned for IGH
- if IGH is zero formally. The ALGOL procedure EXC contained in CBALANCE
- appears in
-
- CCCCBBBBAAAALLLL in line. (Note that the ALGOL roles of identifiers K,L have been
-
-
-
- PPPPaaaaggggeeee 1111
-
-
-
-
-
-
- ____CCCCBBBBAAAALLLL((((3333FFFF)))) ____CCCCBBBBAAAALLLL((((3333FFFF))))
-
-
-
- reversed.) Questions and comments should be directed to B. S. Garbow,
- APPLIED MATHEMATICS DIVISION, ARGONNE NATIONAL LABORATORY
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
-
- PPPPaaaaggggeeee 2222
-
-
-
-